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Abstract 

We study the hadron-quark phase transition in the interior of neutron stars. The relativistic 
mean field (RMF) theory is adopted to describe the hadronic matter phase, while the Nambu-Jona- 
Lasinio (NJL) model is used for the quark matter phase. The influence of the hadronic equation of 
state on the phase transition and neutron star properties are investigated. We find that a neutron 
star possesses a large population of hyperons, but it is not dense enough to possess a pure quark 
core. Whether the mixed phase of hadronic and quark matter exist in the core of neutron stars 
depends on the RMF parameters used. 
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I. INTRODUCTION 



Neutron stars are laboratories for dense matter physics, since they contain matter in 
one of the densest forms found in the universe. It is expected that the deconfmement phase 

n 

transition may occur in the core of massive neutron stars [1[ . The study of the hadron-quark 
phase transition at high density is of great interest in both nuclear physics and astrophysics. 
It has been pointed out by Glendenning [2J] that the hadron-quark phase transition in neutron 
stars may proceed through a mixed phase of hadronic and quark matter over a finite range 
of pressures and densities according to the Gibbs criteria for phase equilibrium. Such phase 
transition has received much attention in neutron star phvsics , [10, [111 • It is 



believed that hyperons may appear around twice normal nuclear matter density through the 
weak interaction [jjj], which usually occur earlier than the hadron-quark phase transition. 
The inclusion of hyperons in the hadronic phase alters the threshold density and properties 
of the mixed phase significantly. In general, the presence of new degrees of freedom, such as 
hyperons and quarks, tends to soften the equation of state (EOS) at high density and lower 
the maximum mass of neutron stars. 

To study the hadron-quark phase transition, we need models to describe hadronic matter 
and quark matter. Unfortunately, there is no single model which can be used to describe both 
phases and the dynamic process of the phase transition. We have to use different approaches 
for the description of the two phases, then perform the Glendenning construction for the 
charge- neutral mixed phase where both hadronic and quark phases coexist [2J . In this work, 
we adopt the relativistic mean field (RMF) theory to describe the hadronic matter phase, 
while the Nambu-Jona-Lasinio (NJL) model is used for the quark matter phase. The choice 
of the NJL model is motivated by the fact that this model can successfully reproduce many 
aspects of quantum chromo dynamics such as the nonperturbative vacuum structure and 
dynamical breaking of chiral symmetry k], 14, ^1. We adopt a three-flavor version of the 
NJL model to describe the quark matter phase 13(. With a definite EOS for quark matter 
based on the NJL model, we examine the influence of the hadronic EOS on the hadron-quark 
phase transition and neutron star properties. 

We use the RMF theory to describe the hadronic matter phase. The RMF theory 



has been success 
clei Il6, 17, 3, U 



ully and widely used for the description of nuclear matter and finite nu- 



201 ] . It has also been applied to provide the EOS of dense matter for the use 
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in supernovae and neutron stars 2l|, |22j . In the RMF approach, baryons interact through the 



exchange of scalar and vector mesons. The meson-nucleon coupling constants are generally 
determined by fitting to some nuclear matter properties or ground-state properties of finite 
nuclei. However, there are large uncertainties in the meson-hyperon couplings due to lim- 
ited experimental data. One can use the coupling constants derived from the quark model, 
or those values constrained by reasonable hyperon potentials. The two additional strange 
mesons, a* and 0, were originally introduced in order to obtain the strong attractive hyperon- 
hyperon (YY) interaction deduced from the earlier measurement [23] . A recent observation 
of the double-A hypernucleus ^He, called the Nagara event 2^], has had a significant im- 
pact on strangeness nuclear physics. The Nagara event provides unambiguous identification 
of ^ A He production with precise AA binding energy value -Baa = 7.25±0.191q u MeV, which 
suggests that the effective AA interaction should be considerably weaker (AI?aa ~ 1 MeV) 
than that deduced from the earlier measurement (AB^ ~ 5 MeV). The weak YY in- 
teraction suggested by the Nagara event has been used to reinvestigate the properties of 
multistrange systems, and it has been found that the chan ge o f YY interactions affects the 



properties of strange hadronic matter dramatically [25|, [26|, l27|]. In order to examine the in- 
fluence of the hadronic EOS on the hadron-quark phase transition, we employ two successful 
parameter sets of the RMF model, NL3 28] and TM1 [29], which have been widely used 
for the description of nuclear matter and finite nuclei including unstable nuclei. For each 
parameter set of the nucleonic sector, we consider two cases of hyperon-hyperon interactions, 
the weak and strong YY interactions. By comparing the results with different parametriza- 
tions in the RMF model, we evaluate how sensitive the hadron-quark phase transition is to 
the hadronic EOS used in the calculation. 

For a comprehensive description of neutron stars, we need not only the EOS at high 
density for the interior region but also the EOS for the inner and outer crusts, where the 
density is low and heavy nuclei exist. For the nonuniform matter at low density, we adopt 



a relativistic EOS based on the RMF theory with a local density approximation 21 



22). 



The nonuniform matter is modelled to be composed of a lattice of spherical nuclei immersed 
in an electron gas with or without free neutrons dripping out of nuclei. As the density 
increases, heavy nuclei dissolve and the optimal state is a uniform matter consisting of 
neutrons, protons, and leptons (electrons and muons) in (3 equilibrium. The low density EOS 
is therefore matched to an EOS of uniform nuclear matter at around 10 14 g/cm 3 [22]. At 
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higher densities, some additional degrees of freedom such as hyperons and quarks may occur, 
and the hadron-quark phase transition could proceed through a mixed phase of hadronic 



uark pnase trar 



and quark matter [5|, |g, 13, |8|, |9|, lid, yj| • Applying the EOS of neutron star matter over a wide 



density range, we study the neutron star properties by solving the Tolman-Oppenheimer- 
Volkoff equation, and examine whether or not the deconfined quark phase can exist in the 
core of neutron stars. 

This paper is arranged as follows. In Sec. [Til we discuss the EOS for hadronic matter in 
the RMF theory, and the parameters used in the calculation. In Sec. Illlj. the NJL model 
is used for the description of quark matter. In Sec. IIVI we investigate the hadron-quark 
phase transition of neutron star matter, and examine the influence of the hadronic EOS. We 
present in Sec. |V]the properties of neutron stars. Section ED is devoted to a summary. 



II. HADRONIC PHASE 



To describe the hadronic matter phase, we adopt the relativistic mean field (RMF) theory, 
in which baryons interact via the exchange of mesons. The baryons considered in this work 
are nucleons (p and n) and hyperons (A, S, and H). The exchanged mesons include isoscalar 
scalar and vector mesons (a and u), isovector vector meson (p), and two additional hidden- 
strangeness mesons (a* and <fi). For neutron star matter consisting of a neutral mixture of 
baryons and leptons in (3 equilibrium, we start from the effective Lagrangian 

C-RMF = 2J -m B - QaBO ~ 9o*B<y* ~ 9wBl^ 

B 

-^BT^f - 9pBl,j,TiPi] ipB + -d^ad^a - ^m 2 a a 2 
-\92cr 3 - \g-^ - \w, v W^ + \mlu^ 

+\d ll o*Po> - \ml,a* 2 - + \m\^^ 

I 

where the sum on B runs over the baryon octet (p, n, A, S + , E°, H°, and the 
sum on I is over electrons and muons (e~ and fi~). In the RMF model, the meson fields are 
treated as classical fields, and the field operators are replaced by their expectation values. 
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The meson field equations in uniform matter have the following form: 
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where a = (a) , u = (u°) , p = (p 30 ) , a* = (a*) , and = (0°) are the nonvanishing 
expectation values of meson fields in neutron star matter. m* B = uib + g a B<J + ga*BV* is the 
effective mass of the baryon species B, and kp is the Fermi momentum. 

In this work, we employ two successful parameter sets of the RMF model, NL3 and TM1, 
as listed in Table [H These parameters have been determined by fitting to some ground-state 
properties of finite nuclei, and they c an p rovide good description of nuclear matter and finite 
nuclei including unstable nuclei 28|, |29j]. As for the meson-hyperon couplings, we take the 
naive quark model values for the vector coupling constants, 



1 
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The scalar coupling constants are chosen to give reasonable hyperon potentials. The poten- 
tial depth of the hyperon species i in matter of the baryon species 7 is denoted by and we 



use U 



(TV) 



-28 MeV, U. 



(TV) 



+30 MeV, and ui N) 



-18 MeV 



30 



31 



321 ] to determine the 



scalar coupling constants g a \, g a -£, and g a s, respectively. The hyperon couplings to strange 
meson a* are restricted by the relation ~ l/j^ ~ 2U^ ~ 2U^ obtained in Ref. 33]. 
We consider two cases of hyperon-hyperon (YY) interactions. The strong YY interaction 



deduced from the earlier measurement 



23|] suggests Up ~ -20 MeV^while the weak YY 
interaction implied by the Nagara event suggests ~ —5 MeV 25, 26, 27]. In Table [III 
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we present the meson-hyperon couplings determined by these hyperon potentials. The hy- 
peron masses are taken to be m A = 1115.7 MeV, ms = 1193.1 MeV, and m s = 1318.1 MeV, 
while the strange meson masses are m a * = 975 MeV and = 1020 MeV. 

For neutron star matter consisting of a neutral mixture of baryons and leptons, the f3 
equilibrium conditions without trapped neutrinos are given by 

Hp = HX+ = Hn ~ He, (8) 

Ha = Hs° = fj, E o = fi n , (9) 
Ht.- = Hz- = Hn + He, (10) 

Hp- — He, (11) 

where /!« is the chemical potential of species i. At zero temperature the chemical potentials 
of baryons and leptons are expressed by 

Hb = yj kf' + rrig + g^B^ + 9^b4> + 9 p bt 3 bP, (12) 
IH = ^k l F 2 + ml (13) 

The charge neutrality condition is given by 

n p + n s + = n e + n M + n s - + n s - , (14) 

where n« = (k l F ) 3 /(3tc 2 ) is the number density of species i. We can solve the coupled 
equations self-consistently at a given baryon density % = n p + n n + + ns+ + n^o + n^- + 
n s o + n s - . The total energy density and pressure of neutron star matter are written by 

1 I 1 11 

£hp = Yl ~2 J \A 2 + m B k2dk + 2 m l a2 + + - A 9^ 

B ^ 

,1 22,3 4,1 22,1 2 *2,1 2/2 

+ 2 m ^ u + 4 C 3^ + 2 m P P + 2 m a* a + 



+ / Jk 2 + m 2 k 2 dk, (15) 



k ? k'dk 1 2 2 1 3 1 4 



,1 22,1 4,1 22 1 2 *2,1 2/2 

+ 2 m ^ +^c 3 uj +-m p p --m^a + -m^ 
1^ 1 k^dk 

+ 3^/ TFT^f- (16) 
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III. QUARK PHASE 



In this section, we adopt a three-flavor version of the N JL model to describe the deconfined 
quark phase. The Lagrangian is given by 

8 

C-njl = q {yy^ - m°) q + [(v^aq) 2 + [qiihKqf] 

a=0 

-K {det [q (1 + 75) q] + det [q (1 - 75) q]} , (17) 

where q denotes a quark field with three flavors (Nf = 3) and three colors (N c = 3). 
m° =diag(m°, m? d , 771°) is the current quark mass matrix, and we assume isospin symmetry 
m° u = m d = m®. The coupling constants G and K have dimension energy -2 and energy -5 , 
respectively. The model has five parameters, namely, the current quark masses m° q and m° s , 
the coupling constants K and G, and the momentum cutoff A. In the present calculation, we 



employ the parameters given in Ref. [34|, m° q = 5.5 MeV, m° s = 140.7 MeV, A = 602.3 MeV, 
GA 2 = 1.835, and KA 5 = 12.36. These parameters have been determined by fitting f w , m n , 
w,k, and m > to their empirical values, while the mass of the 77- meson is underestimated by 
about 6% |ia|. 

In the NJL model, the quark gets constituent quark mass by spontaneous chiral symmetry 
breaking. The constituent quark mass in vacuum is much larger than the current quark 
mass mj. In the quark matter at high density, the constituent quark mass m* becomes 
approximately the same as m", which reflects the restoration of chiral symmetry. Within 
the mean-field approximation, m* is obtained by solving the gap equation 



m\ - 4G(q iqi ) + 2K{q jqj ){q k q k ), (18) 



with (i,j,k) being any permutation of (u,d,s). The quark condensate Cj = (g^i) is given 
by 



3 f A m* 



d = / ; I k 2 dk, (19) 

where k l F denotes the Fermi momentum of the quark flavor i, which is connected with the 
number density n» and the chemical potential /!« via 

\3 



n t = (20) 



IH = dk^ + mf. (21) 
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The energy density of the quark system is given by 



£njl 



E 



3 ' A 



7T 2 



I ^Jk 2 + mf k 2 dk 



+2G {C 2 U + C\ + Cl) - AKC u C d C s - e , (22) 

where Eq is introduced to ensure Enjl = in the vacuum. 

For the quark matter consisting of a neutral mixture of quarks (u, d, and s) and leptons 
(e and //) in (5 equilibrium, the charge neutrality condition is expressed as 

2 1 

-n u ~ - (n d + n s ) -n e -n^ = 0, (23) 
the j3 equilibrium conditions are given by 

A* s = fJ>d = fJ-u + He, (24) 
Hp = ii e - (25) 

The coupled equations can be solved self-consistently at a given baryon density % = 
(n u + rid + n s ) /3. The total energy density and pressure including the contributions from 
both quarks and leptons are given by 

sqp = E NJL + V \ [ F Jk 2 + mf k 2 dk, (26) 

P QP = E n iVi~ £ QP- (27) 

IV. HADRON-QUARK PHASE TRANSITION 

In this section, we study the hadron-quark phase transition which may occur in the core 
of massive neutron stars. It has been discussed extensively in the literature that a mixed 
phase of hadronic and quark matter could exist over a finite range of pressures and densities 
according to the Gibbs criteria for phase equilibrium. In the mixed phase, the local charge 
neutrality condition is replaced by a global one. This means that both hadronic and quark 
matter are allowed to be separately charged. The condition of global charge neutrality is 
expressed as 

XnT + (1 - X) nf p = 0, (28) 
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where x is the volume fraction occupied by quark matter in the mixed phase, which increases 
from x — in the pure hadronic phase to x = 1 i n the pure quark phase. n^ p and nf p 
denote the charge densities of hadronic phase and quark phase, respectively. Without the 
constraint of local charge neutrality, we impose that the two phases are in weak equilibrium 
and described by two independent chemical potentials (/x n , /i e ). The Gibbs condition for 
phase equilibrium at zero temperature is then given by 

PHP {Hn, Me) = PQP (Mm Me) • (29) 

Using Eq. ([29"]) we can calculate the equilibrium chemical potentials of the mixed phase 
where Pup = Pqp = Pmp holds. The energy density and the baryon density in the mixed 
phase are given by 

£mp = X£qp + (1 - X) £ffp, (30) 

and 

nr = Xn Q B P + (l- X )nr- (31) 

We show in Fig. [1] the possible phase structure of neutron star matter using the RMF 
model for the hadronic phase and the NJL model for the quark phase. To examine the influ- 
ence of the hadronic EOS on the hadron-quark phase transition, we employ two successful 



parameter sets of the RMF model, NL3 28j and TM1 29|, with both the weak and strong 
YY interactions. The shaded regions correspond to the mixed phase. It is shown that a pure 
hadronic phase is favored at low density. The mixed phase appears at the critical density 
n B where the pressure of the pure hadronic phase becomes to be lower than the pressure 
of the mixed phase. The fraction of quark matter x increases with increasing density in the 

(2) 

mixed phase. It turns to be a pure quark phase at the critical density n B where the pressure 
of the pure quark phase is above the pressure of the mixed phase. The critical densities n B 
and n B are sensitive to the RMF parameters used, and we get n B ' ~ 0.50 fm" d (0.71 fm" d ) 
and rS B ~ 0.75 fm -3 (0.98 fm~ 3 ) for the NL3 parameter set with the weak (strong) YY 
interaction, while n B ~ 1.31 fm~ (1.75 fm~ ) and n B ~ 2.03 fm~ (2.49 fm - ) are ob- 
tained for the TM1 cases. It is found that the hadron-quark phase transition occurs at lower 
densities in the NL3 model than in the TM1 model, and in general the weak YY interaction 
leads to an earlier appearance of the mixed phase. In order to estimate the influence of 
hadronic EOS on the deconfinement phase transition, we plot in Fig. [2] the hadronic EOS 
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in the different cases and the quark EOS in the NJL model with local charge neutrality as 
a function of the neutron chemical potential \x n . The crossing of the hadronic EOS with 
the quark EOS marks the transition point between the charge neutral hadronic matter and 
quark matter. It is seen that the NL3 model favors the phase transition at a lower \i n . A 
more realistic treatment of the phase transition is to release the constraint of local charge 
neutrality, which leads to the existence of the mixed phase of charged hadronic and quark 
matter over a finite range of pressures and densities as shown in Fig. [TJ In general, a harder 
hadronic EOS also favors an earlier appearance of the mixed phase. 

In Fig. [3] we plot the full EOS in the form P = P(e), which consists of three parts: (a) 
the charge neutral hadronic matter phase at low density (b) the mixed phase of charged 
hadronic and quark matter (c) the charge neutral quark matter phase at high density. The 
mixed phase part of the EOS is shaded gray, where the pressure varies continuously. It is 
shown that the onset and width of the mixed phase depend on the RMF parameters used 
in the calculation. The NL3 model leads to earlier appearance of the mixed phase than the 
TM1 model, and the weak YY interaction favors earlier onset of the mixed phase than the 
strong YY interaction. This is mainly because that a harder hadronic EOS prefers an earlier 
hadron-quark phase transition. By comparing the results of different cases, we can see the 
influence of the hadronic EOS on the hadron-quark phase transition. 

We present in Figs. H] and [5] the particle fraction = n^/Ve as a function of the total 
baryon density tib- At low densities the fractions Y p , Y e , and Y" M increase with increasing 
density. When the Fermi energy of nucleons exceeds the rest mass of hyperons, the conversion 
of nucleons to hyperons is energetically favorable, and it can relieve the Fermi pressure 
of nucleons. The fraction of hyperons increases with increasing density before the mixed 
phase occurs. Quarks appear at the critical density , then the fractions Y u , Yd, and Y s 
increase rapidly with increasing density. The hadronic matter completely disappears at 

(2) 

the critical density n B where the pure quark phase occurs. At extremely high density, 
Y u ~Y^~Y^~l/3 due to the restoration of chiral symmetry. It is shown that the 
composition of neutron star matter depends on the RMF parameters used in the calculation. 
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V. NEUTRON STARS PROPERTIES 



In this section, we investigate the properties of neutron stars by solving the Tolman- 
Oppenheimer-Volkoff equation with the EOS over a wide density range. For the nonuniform 
matter at low density, which exists in the inner and outer crusts of neutron stars, we adopt a 



relativistic EOS based on the RMF theory with a local density approximation [2l|, |22|]. The 
nonuniform matter is modelled to be composed of a lattice of spherical nuclei immersed in 
an electron gas with or without free neutrons dripping out of nuclei. The low density EOS 
is matched to the EOS of uniform hadronic matter at the density where they have equal 
pressures. The pure hadronic phase ends at the critical density rig , and the pure quark 

(2) 

phase starts at the critical density n B . The values of these critical densities depend on 
the RMF parameters used. The neutron star properties are mainly determined by the EOS 
at high density. We calculate neutron star profiles in order to examine whether or not the 
deconfined quark phase can exist in the core of neutron stars. 

In Fig. [6] we present the mass-radius relation using the EOS of the NL3 and TM1 models 
with both the weak and strong YY interactions. It is shown that the results depend on the 
RMF parameters used in the calculation. Since the pressure and density inside neutron stars 
decrease from the center to the surface, the most possible region where the deconfined quark 
phase can exist is the center of the neutron star with maximum mass. We list in Table III 1 1 
the properties of neutron stars with the maximum mass. It is found that the central baryon 
density n c is between and for the NL3 model. This means that the neutron star can 
possess a mixed phase core, but it is not dense enough to possess a pure quark core. On the 
other hand, the values of n c in the TM1 model are smaller than rig, which means that the 
neutron star is only composed of hadronic matter. As can be seen in Table IHfl the neutron 
star properties significantly depend on the RMF parameters used. 



VI. CONCLUSIONS 



We have studied the hadron-quark phase transition at high density, which may occur 
in the core of massive neutron stars. In the present work, we have adopted the RMF 
theory to describe the hadronic matter phase, while the NJL model has been used for the 
quark matter phase. With a definite EOS for the quark phase, we examine the influence of 
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the hadronic EOS on the deconfinement phase transition and neutron star properties. In 
this paper, we employ two successful parameter sets of the RMF model, NL3 and TM1, 
which have been widely used for the description of nuclear matter and finite nuclei including 
unstable nuclei. For each parameter set of the nucleonic sector, we consider two cases of 
hyperon-hyperon interactions, the weak and strong YY interactions. The hadron-quark 
phase transition can proceed through a mixed phase of hadronic and quark matter over a 
finite range of pressures and densities according to the Gibbs criteria for phase equilibrium. 
We have found that the mixed phase starts at rig ~ 0.50 fm -3 (0.71 fm -3 ) and ends 
at n B ~ 0.75 fm - (0.98 fm~ ) for the NL3 parameter set with the weak (strong) YY 
interaction, while rig ~ 1.31 fm~ 3 (1.75 fm~ 3 ) and rig ~ 2.03 fm~ 3 (2.49 fm -3 ) have been 
obtained for the TM1 cases. It is shown that the hadron-quark phase transition occurs 
at lower densities in the NL3 model than in the TM1 model, and in general the weak YY 
interaction leads to an earlier appearance of the mixed phase. By comparing the results with 
different parametrizations in the RMF model, we can see how sensitive the deconfinement 
phase transition is to the hadronic EOS used in the calculation. 

We have calculated the properties of neutron stars using the EOS over a wide density 
range. The star properties such as their masses and radii are mainly determined by the 
EOS at high density. We have found the maximum mass of neutron stars falls in the range 
1.62 ~ 2.05 M Q for the RMF parameters used. With the NL3 model, the mixed phase can 
exist in the core of massive neutron stars, but no pure quark phase can exist. For the TM1 
model, the neutron star is not dense enough to possess the mixed phase, and therefore the 
hadron-quark phase transition could not occur inside neutron stars in this approximation. 
We conclude that the maximum mass and composition of neutron stars depend on the 
hadronic EOS adopted in the calculation. 

It is very interesting to compare our results with those previously published in the litera- 
ture the authors studied the possible hadron-quark phase transition 
utilizing the same NJL model as used in the present work to describe the deconfined quark 
phase, while the hadronic phase was described by several RMF models (TM1, TM2, GL85, 
and GPS as listed in Table I of Ref. p]). Comparing with our cases, different hyperon 
potentials were used to constrain the hyperon couplings, but some of their values are not 
supported by recent experimental observations. They found that the use of the GPS model 
for the hadronic phase leads to the onset of the mixed phase at the energy density e w 7eo 
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(eo = 140 MeV/fm 3 is the normal nuclear energy density), while the mixed phase does not 
appear below e « 10eo for the TM1, TM2, and GL85 models. They concluded that within 
the model constructed in their calculation the appearance of deconfined quark matter in 
the center of neutron stars turns out to be very unlikely. This is consistent with our results 
although the hyperon couplings are different between our study and theirs. In Ref. [6], the 
authors employed an extended MIT bag model to describe the quark phase and four RMF 
models (TM1, TM2, GL85, and GPS) to describe the hadronic phase. They studied the 
influence of different hadronic EOS and the influence of the model parameters of the quark 
phase on the properties of the phase transition. Using the MIT bag model, they found the 
mixed phase could appear below e ~ 2eo, which is much lower than the threshold density 
obtained using the NJL model. This can be seen by comparing Fig. 4 of Ref. with Fig. 9 
of Ref. 5} and Fig. 3 of the present paper. In Ref. [?], the authors investigated the de- 
confinement phase transition at zero and finite temperature using both the MIT bag model 
and the NJL model for describing the quark phase and a RMF model (GL model) for the 
hadronic phase. For the NJL model, however, they adopted another parameter set which 
favors an earlier onset of the mixed phase than the parameter set used in our calculation. 
For the hyperon couplings, they used three choices and verified that for the NJL model the 
onset of the mixed phase is very sensitive to the choice of the hyperon couplings. Using the 
NJL model for the quark phase, they obtained two different behaviors. One is the hyperons 
appear before the quarks and the mixed phase occurs at much higher densities (~ 5po with 
the normal nuclear matter density po — 0.153 fm~ 3 ). The other is the quarks appear at 
lower densities than the hyperons and the mixed phase occurs at ~ 2po- For almost all 
EOS used in their study, the central energy density of the maximum-mass neutron star falls 
inside the mixed phase, so the star can contain a core constituted by a mixed phase, but 
it is not dense enough to possess a pure quark core. In Ref. [8j, the authors investigated 
the deconfinement phase transition using an effective-field-theory-motivated RMF (E-RMF) 
model for the hadronic phase and considering both unpaired quark matter (UQM) described 
by the MIT bag model and paired quarks described by the color-flavor locked (CFL) phase 
for the quark phase. They mentioned that they could not get any mixed phase with the 
original G2 parameter set in the E-RMF model, so they changed incompressibility K from 
215 to 300 MeV and the effective mass m* N /m^ from 0.664 to 0.7 to determine a modified 
parameter set G2* as given in Table I of Ref. [8J . For the hyperon couplings, they assumed 
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that all the hyperons in the octet have the same couplings. Using the G2* parameter set in 
the E-RMF model for the hadronic phase and B 1 ^ = 170 MeV in the UQM model for the 
quark phase, they obtained the onset of the mixed phase at ~ l-3p 0) while the appearance 
of the mixed phase starts at ~ 2.4p using the CFL model with B 1 ^ = 188 MeV for the 
quark phase. In our cases, the mixed phase occurs at ~ 3.3p (8.6p ) using the NL3 (TM1) 
model with the weak YY interaction for the hadronic phase and the NJL model for the 
quark phase. For simplicity, we have not considered the color-flavor locked phase in the 
present work. By comparing all these results, we conclude that the deconfmement phase 
transition is very sensitive to both the hadronic EOS and the quark EOS, whether or not 
the deconfmed quark matter appears in the center of neutron stars depends on the models 
adopted in the calculation. 
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TABLE I: The parameter sets NL3 [28| and TM1 



291 ] used in the calculation. The masses are 



given in MeV. 

m N m a m p g aN g^ N g pN g 2 (fm _1 ) g 3 c 3 

NL3 939.0 508.194 782.501 763.0 10.217 12.868 4.474 -10.431 -28.885 
TM1 938.0 511.198 783.0 770.0 10.029 12.614 4.632 -7.233 0.618 71.308 



TABLE II: The scalar coupling constants for the hyperons. 

9a A QctY, 9aE 9a* A(S) 9v*E 

NL3 (weak YY) 6.269 4.709 3.242 5.595 11.765 

NL3 (strong YY) 6.269 4.709 3.242 7.138 12.809 

TM1 (weak YY) 6.170 4.472 3.202 5.412 11.516 

TM1 (strong YY) 6.170 4.472 3.202 7.018 12.600 



TABLE III: The properties of neutron stars with the maximum mass M max . The central energy 
density, pressure, and baryon number density are denoted by e c , P c , and n c , respectively. R and 
-Rmp denote the radii of the star and its mixed phase core. 

e c P c n c R Rmp 

(M ) (10 15 g/cm 3 ) (10 35 dyn/cm 2 ) (fm~ 3 ) (km) (km) 

NL3 (weak YY) 2.02 1.42 2.08 0.69 13.80 4.21 

NL3 (strong YY) 2.05 1.46 2.62 0.71 13.48 0.31 

TM1 (weak YY) 1.71 1.50 2.08 0.75 13.07 - 

TM1 (strong YY) 1.62 1.23 1.32 0.63 13.62 - 
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FIG. 1: (Color online) The pressure P as a function of the baryon density The shaded regions 
correspond to the mixed phase (MP). The dashed and dot-dashed lines show the pressures of 
hadronic phase (HP) and quark phase (QP), respectively. 
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FIG. 2: (Color online) The pressure P as a function of the neutron chemical potential fj, n for 
hadronic and quark matter with local charge neutrality. The crossing of the hadronic EOS with 
the quark EOS (NJL) marks the transition point between the charge neutral hadronic matter and 
quark matter. 
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FIG. 3: (Color online) The full EOS of neutron star matter in the form of pressure P versus energy 
density e. The shaded regions correspond to the mixed phase (MP). The dashed and dot-dashed 
lines show the pressures of hadronic phase (HP) and quark phase (QP), respectively. 




FIG. 4: (Color online) The particle fraction Yi = rii/nB as a function of the total baryon density 
riB for the NL3 model. 
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FIG. 5: (Color online) The particle fraction Y( = m/riB as a function of the total baryon density 
ub for the TM1 model. 




FIG. 6: (Color online) The mass-radius relation for neutron stars. 
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